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Abstract —In this paper, we study the resource allocation 
algorithm design for energy-efficient simultaneous wireless in¬ 
formation and power transfer (SWIPT) systems. The considered 
system comprises a transmitter, an information receiver, and 
multiple energy harvesting receivers equipped with multiple 
antennas. We propose a multi-objective optimization framework 
to study the trade-off between the maximization of the energy 
efficiency of information transmission and the maximization 
of wireless power transfer efficiency. The proposed problem 
formulation takes into account the per antenna circuit power 
consumption of the transmitter and the Imperfect channel state 
information of the energy harvesting receivers. The adopted non- 
convex multi-objective optimization problem is transformed into 
an equivalent rank-constrained semidefinite program (SDP) and 
optimally solved by SDP relaxation. Numerical results unveil an 
interesting trade-off between the considered conflicting system 
design objectives and reveal the benefits of multiple transmit 
antennas for improving system energy efficiency. 

I. Introduction 

In recent years, the development of wireless communication 
networks worldwide has triggered an exponential growth in 
the number of wireless devices and sensors for applications 
such as e-health and environmental monitoring. The related 
tremendous increase in the number of transmitter(s) and re- 
ceiver(s) has also led to a huge demand for energy and a better 
energy management. Hence, energy efficient system designs, 
which adopt energy efficiency (bit-per-Joule) as the perfor¬ 
mance metric, have been recently proposed Q-ll. In d, 
energy-efficient power allocation schemes were proposed for 
cognitive radio systems. In d, energy-efficient link adaptation 
was investigated for the maximization of energy efficiency 
in frequency-selective channels. In d, the authors proposed 
a resource allocation algorithm design for energy-efficient 
communication in multicarrier communication systems with 
hybrid energy harvesting base stations. Although energy- 
efficient resource allocation algorithm designs for traditional 
communication networks have been studied in the literature, 
mobile receivers are often powered by batteries with limited 
energy storage which remain the system performance bottle¬ 
necks in perpetuating the lifetime of wireless networks. 

Energy harvesting (EH) based communication system de¬ 
sign is a viable solution for prolonging the lifetime of energy- 
limited devices. Conventional natural sources, such as wind, 
solar, and biomass, have been exploited as energy sources 
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for fixed-location, outdoor transmitters. However, these natural 
energy sources are often location and weather dependent and 
may not be suitable for mobile receivers. On the other hand, 
wireless power transfer (WPT) via electromagnetic waves in 
radio frequency (RE) enables a comparatively controllable 
energy harvesting for mobile receivers. In fact, recent progress 
in the development of RE-EH circuitries has made RE-EH 
practical for low-power consumption devices a-m, e.g. 
wireless sensors. Besides, WPT enables the dual use of the 
information carrier for simultaneous wireless information and 
power transfer (SWIPT) lISl- lfTOl . Different from traditional 
wireless communication systems, where data rate and energy 
efficiency are the most fundamental system performance met¬ 
rics, in SWIPT systems, the wireless energy transfer efficiency 
is an equally important QoS metric. Thus, the design of 
resource allocation algorithms should take into account the 
emerging need for energy transfer efficiency. In ifTTl . the 
authors studied the fundamental rate-energy trade-off region 
for optimal information beamforming. In lIT^ . power alloca¬ 
tion, user scheduling, and subcarrier allocation were jointly 
designed to enable an energy-efficient multicarrier SWIPT 
system. In ca, the authors proposed the use of large scale 
multiple-antenna systems for improving energy efficiency of 
SWIPT. Although energy-efficient data communication de¬ 
sign and energy-efficient WPT have already been studied 
individually, the trade-off between these two system design 
paradigms is still unclear for SWIPT systems. In particular, 
these two design goals may conflict with each other but both 
are desirable to system designer. However, the single-objective 
resource allocation algorithms proposed in E-ii, m-ini, 
lfT3]l may no longer be applicable in energy-efficient SWIPT 
networks. 

In this paper, we address the above issues. To this end, 
we formulate the resource allocation algorithm design as a 
multi-objective optimization problem which strikes a balance 
between the maximization of energy efficiency of information 
transmission and the maximization of WPT efficiency. The 
resulting non-convex optimization problem is solved optimally 
by semidefinite programming (SDP) relaxation. Simulation 
results illustrate the trade-off between the conflicting system 
design objectives. 

II. System Model 

In this section, we first define the adopted notations and 
then present the channel model for energy-efficient SWIPT 



Fig. 1. A downlink communication system with an information receiver (IR) 
and J = 2 energy harvesting receivers (ERs). 

networks. 

A. Notation 

, Tr(A), and Rank(A) represent the Hermitian trans¬ 
pose, trace, and rank of matrix A; A ^ 0 indicates that A is 
a positive semidefinite matrix; matrix I^v denotes an N x N 
identity matrix. vec(A) denotes the vectorization of matrix A 
by stacking its columns from left to right to form a column 
vector. A (g) B denotes the Kronecker product of matrices A 
and B. [B]a:b_c:d returns the a-th to the 6-th rows and the c-th 
to the d-th columns block submatrix of B. denotes 

the space of N x M matrices with complex entries. 
represents the set of all N-hy-N complex Hermitian matrices. 
The distribution of a circularly symmetric complex Gaussian 
(CSCG) vector with mean vector x and covariance matrix 
S is denoted by CAf{K, S), and ~ means “distributed as”. 
£{■} denotes statistical expectation. H-H and H-Hi? denote the 
Euclidean norm and the Frobenius norm of a vector/matrix, 
respectively. Re(-) extracts the real part of a complex-valued 
input. 

B. Channel Model 

We focus on a downlink SWIPT system. The system con¬ 
sists of a transmitter, a single-antenna information receiver 
(IR), and multiple energy harvesting receivers (ERs). The 
transmitter is equipped with At antennas and each ER is 
equipped with Ar receiving antennas. We assume that the 
ERs are roaming wireless terminals from other communication 
systems searching for additional power supply in the RE. The 
transmission is divided into time slots. In each time slot, the 
transmitter sends a precoded information signal and an energy 
signal simultaneously to facilitate information transmission to 
the IR and power transfer to the ERs, cf. Figure[T] The transmit 
signal is given by 

X = WlS-|-WE, (1) 

where s € C is the information-bearing symbol with 
a{|s|^} = 1 and wi S is the corresponding informa¬ 
tion beamforming vector. We G is the energy signal 

facilitating energy transfer to the ERs. The energy signal we 
is a deterministic pseudo-random sequence with zero mean 
and covariance matrix We. Since we is generated at the 
transmitter by a pseudo-random sequence generator with a 
predefined seed, the energy signal is known to the IR. Thus, 


the interference caused by the energy signal can be completely 
cancelled at the IR. 

We assume a narrow-band slow fading channel between the 
transmitter and receivers. Then, the received signals at the IR 
and ER j are expressed as 

yiR ^ h^(wis-f We)- f ni, (2) 

y™ = Gf (wis-f We)- f n®, (3) 

where h G is the channel vector between the trans¬ 
mitter and the IR, and Gj € is the channel matrix 

between the transmitter and ER j. Variables h and Gj capture 
the joint effect of multipath fading and path loss. ni G C and 
e C^RXi are additive white Gaussian noises (AWGNs) at 
the IR and ER j, respectively, and are distributed as CA(0, cr^) 
and CA(0, ctIItvr). 

III. Resource Allocation Algorithm Design 

In this section, we present the adopted performance metrics 
and the problem formulation. 


A. Achievable Rate, Harvested Energy, and Energy Efficiency 
The achievable rate (bit-per-second) at the IR is given by 
/ W^^Hwt\ 

i? = Blog2(n- ^ -), (4) 

where B is the system bandwidth and H = hh^. We 
note that the interference caused by the energy signal, i.e., 
Tr(wgHwE), is removed from the IR via successive in¬ 
terference cancellation before the IR decodes the desired 
information since the energy signal we is known to the 
receiver. On the other hand, both the information signal and 
the energy signal can act as RE energy sources for the ER due 
to the broadcast nature of wireless channels. As a result, the 
total harvested energ}{^ at ER j is given by 

ph-v = p,Tr(Gf(wiwf+W e)G,), ( 5 ) 


and r/j is the energy conversion efficiency of ER j which is 
a constant with 0 < pj < 1 and models the energy loss of 
the process of converting the received RE energy to electrical 
energy for storage. We ignore the thermal noise at the receiving 
antenna for energy harvesting as it is relatively small compared 
to the received signal power. 

Energy efficiency is a fundamental system performance met¬ 
ric in modern communication networks. To design a resource 
allocation algorithm for energy-efficient communication, the 
total power consumption has to be included in the optimization 
objective function. Thus, we model the power dissipation 
(Joule-per-second) in the system as 


|!wi|i2+Tr(AVE) , ^ 

-ftot — -- ^ ^ 

where Pb = At Pant + Pc- 


(6) 

(7) 


0 < ^ < 1 is the constant power amplifier efficiency. The 
first term in (j^ is the total power consumption in the power 

’We note that a normalized energy unit, i.e., Joule-per-second, is adopted. 
Therefore, the terms “power” and “energy” are used interchangeably in this 
paper. 















amplifier. N^Pant in 0 accounts for the dynamic circuit 
power consumption which is proportional to the number of 
transmit antennas. Pant denotes the power dissipation at each 
transmit antenna, including the dissipation in the transmit 
filter, mixer, frequency synthesizer, digital-to-analog converter 
(DAC), etc. Pc denotes the fixed circuit power consumption 
for baseband signal processing. 

Therefore, the achievable rate energy efficiency (AR-EE) 
and the energy transfer energy efficiency (ET-EE) of the 
considered system are defined as 

Pl0g2(l+^^^) 

(llwiP + Tr(WE))/^ + iVxPant + Pc 

^harv 

(llwiP + Tr(WE))/^ + iVxPant + Pc ’ 
respectively, where is given in 0. 


‘fiR = 

‘&EH = 


B. Channel State Information (CSI) 

In this paper, we focus on a Time Division Duplex (TDD) 
communication system with slowly time-varying channels. At 
the beginning of each time slot, handshaking is performed 
between the transmitter and the IR. As a result, the downlink 
CSI of the IR can be obtained by measuring the uplink training 
sequences embedded in the handshaking signals. Thus, we 
assume that the transmitter-to-IR fading gain, h, can be 
estimated perfectly at the transmitter. On the other hand, the 
ERs may not directly interact with the transmitter Besides, the 
ERs may be silent for long periods of time. As a result, the CSI 
of the ERs can be obtained only occasionally at the transmitter 
when the ERs communicate with the transmitter. Hence, 
the CSI for the ERs may be outdated when the transmitter 
performs resource allocation. We adopt a deterministic model 
oa, ina to characterize the impact of the CSI imperfection 
for resource allocation design. The CSI of the link between 
the transmitter and ER j is modeled as 


Gj = Gj+AGj, Vj e J}, and (10) 


where Gj G c^txa^r matrix CSI estimate of the chan¬ 

nel of ER j that is available at the transmitter. AG_, represents 
the unknown channel uncertainty and the continuous set 
in (111 defines a continuous space spanned by all possible 
channel uncertainties. Constant £j represents the maximum 
value of the norm of the CSI estimation error matrix AG, for 
ER J. 


IV. Problem Eormulation and Solution 
A. Problem Formulation 

In SWIPT systems, AR-EE maximization and ET-EE max¬ 
imization are both desirable system design objectives. In this 
section, we first propose two problem formulations for single¬ 
objective system design for SWIPT. Each single-objective 
problem describes one important aspect of the system design. 
Then, we consider both system design objectives jointly via 
the multi-objective problem formulation. 


The first system design objective is the maximization of 
AR-EE without the consideration of energy harvesting. The 
corresponding optimization problem is formulatecQ as 


Problem 1: 

AR-EE Maximization: 

maximize 

$IR (12) 

WeGH”t,wi 

subject to 

Cl: ||wi||2-fTr(WE)<Pn,ax, 


C2 : We A 0 


^max in constraint Cl denotes the maximum transmit power 
budget. In addition, covariance matrix We is a positive 
semidefinite Hermitian matrix as indicated by constraint C2. 

The second system design objective is the maximization of 
the ET-EE. The corresponding problem formulation is given 
as 


Problem 2: ET-EE Maximization: 


maximize 

min $Eu 

(13) 

WbGH”t,wi 

AG, 64', 

subject to 

Cl, C2. 



Eor the sake of notational simplicity, we denote the ob¬ 
jective functions in the above problems as Fn, n = 1,2. In 
practice, these two system design objectives are both desirable 
from the system operator perspective. However, it is expected 
that there is a non-trivial trade-off between these objectives. 
In order to meet these conflicting system design objectives 
systematically and simultaneously, we adopt the weighted 
Tchebycheff method for the multi-objective optimization m 
which can provide the complete Pareto optimal set by varying 
predefined preference parameters. To this end, we incorporate 
the two individual system design objectives into a multi¬ 
objective optimization problem (MOOP), which is formulated 
as 


Problem 3: Multi-Objective Optimization Problem: 


minimize 

WeGH*'t,wi 

subject to 


Cl, C2, 


(14) 


where F* is the optimal objective value with respect to 
Problem i. uji is a weight imposed on objective function i 
subject to 0 < Wi < 1 and which indicates 

the preference of the system designer for the i-th objective 
function over the others. In the extreme case, when oJi = 1 and 
uJn = 0,Vn f i. Problem]^ is equivalent to single-objective 
optimization problem i. 


^We note that the considered problem formulation can be easily extended 
to the case with a minimum data rate requirement. Yet, a stringent data rate 
requirement does not facilitate the study the trade-off between different system 
objectives due to the resulting smaller feasible solution set. 









V. Optimization Solutions 
It can be observed that the objective functions of Problems 
m are non-convex functions. In general, there is no well- 
known systematical approach for solving non-convex opti¬ 
mization problems. In order to obtain a tractable solution, 
we first transform the non-convex objective functions using 
the Charnes-Cooper transformation. Then, we use semidefinite 
programming relaxation (SDR) to obtain the resource alloca¬ 
tion solution for the reformulated problem. 

We first reformulate the aforementioned three optimization 
problems by defining a set of new optimization variables; 

Wi = wiwf, e = —, 

-^tot 

Wi = 6iWi, and We = 6»We. (15) 

Then, the original problems can be rewritten with respect to the 
new optimization variables {Wj, We, 0}. Problem 1 becomes 


Problem 4: Transformed AR-EE Maximization Prob¬ 
lem: 


maximize 

Wi,WEGH™T,e 

subject to 


6»log2 (|l + 


Tr(HWi) 


6»cr2 


(16) 


Cl 

C2 

C3 

C4 

C5 


Tr(Wi -f We) < OP, 
Wi ^ 0, We P 0, 
Rank(Wi) < 1, 
Tr(Wi + We) 




+ OP^ < 1, 


0 > 0 , 


where Wi ^ 0,Wi S and C3 are imposed to 

guarantee that Wi = 0wiw^. Constant B is dropped from 
the objective function in Problem since it is independent of 
the optimization variables. Similarly, Problem 2 becomes 


Problem 5: Transformed ET-EE Maximization Prob¬ 
lem: 


_maximize 

Wi,WEGH”T,e,7j. 

subject to 


i=i 

CT-C5. 


C6 : 7 , < min P) 


harv 


(17) 


Vj, 


where 7 j are auxiliary optimization variables. 
Einally, Problem 3 can be written as 


Problem 6: Transformed MOOP: 


minimize 

T 

(18) 

Wi,WEGH”T,e,r 



subject to 

CT-C6, 



C7: uj,iF*-F,)< t, 

i G {4,5}, 


where t is an auxiliary optimization variable. 

Proposition 1: The Problems |4]|^ are equivalent transfor¬ 
mations of the original ProblemsTO^ respectively. 

Proof: The transformation is based on Charnes-Cooper 
transformation. Due to the space limitation, we refer to da 
for proof for a similar problem. ■ 

We note that Problem is a generalization of Problems 4 
and 5. If Problem]^ can be solved optimally by an algorithm, 
then the algorithm can also be used to solve Problems 4 and 
5. Thus, we focus on the method for solving Problem 
It is evident that Problem is non-convex due to the rank- 
one beamforming matrix constraint C3 : Rank(Wi) < 1. 
Besides, constraint C6 involves infinitely many constraints due 
to the continuous uncertainty set d/j . Next, we introduce a 
Lemma which allows us to transform constraint C6 into a 
finite number of linear matrix inequalities (LMIs) constraints. 


Lemma 1 (S-Procedure snu): Let a function fm (x), m G 
{1, 2}, X e be defined as 

/m(x) = x^A^x + 2Re{b^x} -f Cm, (19) 

where G H^, G and Cm G K. Then, the 

implication /i (x) < 0^/2 (x) < 0 holds if and only if there 
exists an w > 0 such that 



Ai 

bi 


A 2 

b2 

OJ 

bf 

Cl 


L^2 

C2_ 


provided that there exists a point x such that /^(x) < 0. 
Now, we apply Lemma 1 to constraint C6. In particular, we 
define gj = vec(Gj), Agj = vec(AGj), Wi = I^Ve d Wi, 
and We = IatrCiWe. By exploiting the fact that |j AG^jlp < 
AgfAgj < e], then we have 

\\AGf\l<e] (21) 

^ C6 : 0 > 7 j-I- min -| Agj^ fWi-f We] Agj 

+2Re{gf (Wi + We) Ag,} + gf (Wi + We) g,}, Vj, 

if and only if there exists a pj > 0 such that the following 
LMIs constraint holds: 

PjTjVtjEWe WEgj _ 

gf We -pj£] - f + gf WEgj 

+ UfWiUg^. A0,Vfc, (22) 

where Ug^ = [ItvrAt ; gj] ■ The new constraint C6 is not only 
an affine function with respect to the optimization variables, 
but also involves only a finite number of constraints. Then, 
we apply the SDP relaxation by removing constraint C3 from 
Problem As a result, the SDP relaxed problem is given by 


C6: = 


^In studying the solution structure of Problem 6, we assume that the optimal 
objective values of Problems 4, 5 are given constants, i.e., F* ,\/p E {4, 5}, 
are known. Once the structure of the optimal resource allocation scheme of 
Problem 6 is obtained, it can be exploited to obtain the optimal solution of 
Problems 4, 5. 














TABLE I 

Simulation Parameters 


Canier center frequency 

915 MHz 

Bandwidth 

200 kHz 

Single antenna power consumption 

Pant = 1 W 

Static circuit power consumption 

P, = 150 W fT^ 

Power amplifier efficiency 

CN 

o 

II 

Transmit antenna gain 

18 dBi 

Noise power 

= —95 dBm 

Transmitter-to-ERs fading distribution 

Rician with Rician factor 6 dB 

Transmitter-to-IR fading distribution 

Rayleigh 

Energy conversion efficiency 

r]j = 0.5 


Problem 7: SDP Relaxed Transformed MOOP: 

minimize 

r (23) 

Wi,WBGH”T,6»,T,7j.pj 


subject to 

Cl, C2, C4, C5,C7, 


C6 : ^ 0,C8 : > 0,Vj, 


which is a convex SDP problem and can be solved by 
numerical convex program solvers such as CVX ifTSll . In 
particular, if the obtained solution Wj of the SDP relaxed 
problem satishes constraint C3, i.e., Rank(Wj) < 1, then it 
is the optimal solution. Now, we study the tightness of the 
SDP relaxation by the following theorem. 

Theorem 1: Assuming that the channels, i.e., h and Gj, are 
statistically independent and Problem|^is feasible, the optimal 
beamforming matrix of Problem is a rank-one matrix with 
probability one, i.e., Rank(Wj) < 1. Besides, for uji > 0, 
the optimal energy signal is = 0. 

Proof: Please refer to the Appendix. ■ 

Therefore, the adopted SDP relaxation is tight. Besides, 
whenever AR-EE is considered, i.e., wi > 0, no dedicated 
energy beam is needed. In fact, the optimal information beam, 
wi, serves as a dual purpose carrier for maximization of 
the energy efficiency of information transmission and WPT 
simultaneously. Eurthermore, Problems 1-2 can be solved by 
SDP relaxation as solving Problem]^ 

VI. Results 

In this section, we present simulation results to demonstrate 
the system performance of multi-objective SWIPT system 
design. The simulation parameters are summarized in Ta¬ 
ble U The IR and J ERs are located 100 meters and 10 
meters from the transmitter. In particular, the ERs are near 
the transmitter with line-of-sight communication channels to 
facilitate energy harvesting. Each ER is equipped with = 2 
antennas for facilitating EH. We assume that the noise powers 
at each antenna of the IR and the ERs are identical, i.e., 
= cr| = tj^. In the sequel, we dehne the normalized 

maximum channel estimation error of ER j as 
with = 0.05, Va, 6 S {1,...,J}. All simulation 

results are obtained by averaging the system performance over 
different multipath channel realizations. 

Eigure depicts the trade-off region for the average ET- 
EE and the average AR-EE achieved by the proposed optimal 



Fig. 2. System performance trade-off region between the average ET-EE 
and the average AR-EE for Pmax = 40 dBm. The double-sided arrows 
indicate the performance gain achieved by the proposed optimal scheme over 
the baseline scheme. 

scheme for different numbers of ERs. The maximum transmit 
power is set to Pmax = 40 dBm. The trade-off region in Eigure 
1^ is obtained by solving Problem 7 via varying the values 
of 0 < Xp < l,Vp S {1,2}, uniformly for a step size of 
0.01 such that J2p = 1- It can be observed that the ET- 
EE is a monotonically decreasing function with respect to the 
AR-EE. The result indicates that AR-EE maximization and 
ET-EE maximization are conflicting system design objectives 
in general. In other words, a resource allocation algorithm 
maximizing the AR-EE cannot maximize the ET-EE simulta¬ 
neously in the considered system. Besides, the trade-off region 
is enlarged for an increasing number of ERs. This is due to 
the fact that a larger portion of the radiated power can be 
harvested when there are more ERs in the system since more 
receivers participate in the energy harvesting process. 

Eor comparison, we also plot the trade-off region of a 
baseline power allocation scheme in Eigure]^ Eor the baseline 
scheme, the covariance matrix of the energy signal We is 
set to zero. Then, maximum ratio transmission (MRT) with 
respect to the IR is adopted for the information beamforming 
matrix Wj. In other words, the beamforming direction of 
matrix Wj is fixed. Then, we optimize the power of Wj 
subject to the constraints in Problem 7. It can be observed that 
the baseline scheme achieves a significantly smaller trade-off 
region compared to the proposed optimal scheme. As a matter 
of fact, the degrees of freedom of the beamforming matrix 
Wi are jointly optimized in our proposed optimal scheme 
via utilizing the CSI of all receivers. On the contrary, the 
information beamformer in the baseline scheme is restricted 
to the range space of the IR. Although the baseline scheme 
is optimal when AR-EE is the only system design objective, 
the information beamformer cannot be steered towards the 
direction of the ERs. Thus, compared to the proposed optimal 
scheme, the baseline scheme is less efficient when ET-EE is 
considered. 

Eigures and illustrate the average AR-EE and the 
average ET-EE versus the total transmit power budget Pmax, 
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Fig. 3. Average AR-EE (bits/joule) versus maximum transmit power budget Eig. 4. Average ET-EE versus maximum transmit power budget Pmax (dBm). 
Pmax (dBm). 


respectively, for J = 10 ERs. The results in Figures and 
are obtained by solving Problem 7 with {oJi = 1 ,uj2 = 0} 
and {wi — 0,uj2 — 1}, respectively. It can be observed from 
Figure that the AR-EE of the proposed optimal scheme 
increases with respect to Pmax monotonically and reaches an 
upper limit where the EE gain due to a higher value of Pmax 
vanishes. This result indicates that once the maximum AR- 
EE is achieved by transmitting a sufficiently large power, any 
additional increase in the transmitted power will incur a loss 
in EE which is prevented by the optimal algorithm. On the 
other hand, it can be seen from Figure that the average 
ET-EE increases slowly in the low transmit power regime but 
increases rapidly in the high transmit power regime. This is 
because for a small transmit power, the ET-EE is dominated 
by the fixed circuit power consumption, Pb, leading to a 
slow increasing rate of ET-EE with respect to the transmit 
power. As the transmit power budget increases, the transmit 
power consumption in the RE becomes significant and the ET- 
EE becomes more sensitive to increases in transmit power 
budget, cf. (|^. On the other hand, the number of transmit 
antennas affects the AR-EE and the ET-EE differently. 
In fact, the maximum AR-EE does not necessarily increase 
with the number of transmit antennas when the per-antenna 
power consumption in considered, cf. Figure]^ This is because 
the AR scales logarithmically with respect to the number of 
transmit antennas. However, the AR gain due to extra transmit 
antennas is not sufficient to compensate the total increased 
energy cost since the circuit power consumption increases 
linearly with respect to N^. Thus, adopting exceedingly large 
numbers of transmit antennas may not be a viable solution for 
information transmission. In contrast, the maximum ET-EE 
increases with Nt as shown in Figure This is due to the 
fact that the ET-EE function in (|^ is a quasi-Iinear function 
with respect to both Wj and We. Thus, a large number of 
transmit antennas is beneficial if ET-EE is the only system 
design objective. 

VII. Conclusions 

In this paper, we studied the resource allocation algorithm 
design for energy-efficient SWIFT networks. The algorithm 


design was formulated as a non-convex MOOP employing 
the weighted Tchebycheff method. The proposed problem 
aimed at the joint maximization of the energy efficiency of 
information transmission and WPT simultaneously. Besides, 
the imperfectness of the CSI of the ERs was also taken into 
account for designing a robust resource allocation algorithm. 
The proposed MOOP was solved optimally by SDP relaxation. 
Simulation results not only reveal the trade-off between the 
studied conflicting systems design objectives, but also shed 
some light on the use of multiple transmit antennas for 
improving system energy efficiency. 


Appendix-Proof of Theorem[T] 

Since Problem [7] satisfies Slater’s constraint qualification 
and is jointly convex with respect to the optimization vari¬ 
ables, strong duality holds. Thus, solving the dual problem is 
equivalent to solving the primal problem. Therefore, Theorem 
[^can be proved by analyzing the dual problem of Problem]^ 
To this end, we define the Lagrangian function C 


= r + a(Tr(Wi+WE)-0P„,ax)-<50 


Vl 


a;i(F*-01og2(l + ^^^^))- 


J 


V2 oj2{F^ - I] 7 j-) - ^ - Tr(XWi) - Tr(YWE) 
i=i 

j J 

Tr(D^^. ) - Y ^^lPJ > (24) 

i=i f=i 


where a, /3, 5, > 0 are dual variables associated with 

constraints Cl, C4, C5, and C8, respectively. Dual variable 
matrices X, Y, and are connected to the EMI constraints 
in C2 and C6, respectively. vi,V 2 are the dual variables for 
constraint C7. 

Then, the dual problem of Problem |7] is given by 


maximize minimize C. (25) 

Q;./3.(5,/ij>0,Y.X,D^ ^0 Wi.WeGH^t, 

^ a — ^ 





































Now, we focus on those Karush-Kuhn-Tucker (KKT) con¬ 
ditions which are useful for the proof: 

Y, X, ^ 0 , a, 13, S, fj.j > 0, Vj, 

ViLOiO 


X = Y- 

Y = (a+^ 
I 

XWi = 0, 
YWe = 0. 


-f Tr(HWi) 

o J Nn ^ 


j=l 1 = 1 




(26) 

H, 


(27) 

F)_ 

a:b,i 

,(28) 

z:d 



(29) 



(30) 


where a = {I — 1 )Nt -1-1,6 = INt, c = {I — 1)Nt + 1, and 
d — INt- Next, we investigate the structure of Wj in the 
following two cases. 

Case 1: For oji = 0 and 102 = 1, AR-EE maximization is 
not considered and ET-EE maximization is the only system 
design objective. Besides, since wi = 0, we have X = Y. 
In other words, Wj has the same functionality as We for 
improving the ET-EE. Thus, without loss of generality and 
optimality, we can set Wi = 0 and Rank(Wi) < 1 holds for 
the optimal solution. 

Case 2: Eor wi > 0, AR-EE maximization is considered 
in the resource allocation algorithm design. Thus, constraint 
C 6 for j = 1 is active, i.e., vi > 0. Besides, from the 
complementary slackness condition in ( |29l l, the columns of 
Wi for the optimal solution lie in the null space of X for 
Wi 7 ^ 0. Therefore, if Rank(X) = — 1, then the optimal 

beamforming matrix must be a rank-one matrix. To reveal the 
structure of X, we show by contradiction that Y is a positive 
definite matrix with probability one. Eor a given set of optimal 
dual variables, {Y, X, Djjg , a,/3, 6 , pj}, and optimal primal 
variables {r, 7 ^, 0, We, Pj}, (25 1 can be written as 


minimize C. 

WiGH”t 


(31) 


Suppose Y is not positive definite, then we can construct 
Wi = rvv^ as one of the optimal solutions of (311, where 


r > 0 is a scaling parameter and v is the eigenvector 
corresponding to one of the non-positive eigenvalues of Y. 
We substitute Wi = rvv^ into ( (3l] l which leads to £ = 
Tr(rYvv^) - r Tr (vv^(X + ,,.X(HW:) )) + ^ 
n is the collection of variables that are independent of Wi. 
Since the channels of Gj and h are assumed to be statistically 
independent, it follows that by setting r — 00 , the dual 
optimal value becomes unbounded from below. Besides, the 
optimal value of the primal problem is non-negative. Thus, 
this leads to a contradiction since strong duality does not hold. 
Therefore, Y is a positive definite matrix with probability one, 
i.e., Rank(Y) = Nt- 

By exploiting ( [27] i and a basic inequality for the rank of 
matrices, we have 


Rank(X) -f Rank 
> Rank (Y) = Nt 


viuiidH 


y0af + Tr(HWi) 


Rank(X) > Nt — Rank ( 


(32) 

(33) 

'’■"‘"*^'1.(34) 


V0(T2-f Tr(HWi)A 


Eurthermore, Wi 7 ^ 0 is required to maximize the energy 
efficiency of data communication. Hence, Rank(X) = A^t — 1 
and Rank(Wi) = 1. Besides, in this case, since Y is full 
rank. We = 0 according to ( [30l l. In other words, utilizing 
only information beam wi is optimal when uji > 0 . 
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